Random local strain effects in homovalent-substituted relaxor ferroelectrics: 
a first-principles study of BaTio.74Zro.26C)3 
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We present first-principles supercell calculations on BaTio.74Zro.26O3, a prototype material for 
relaxors with a homovalent substitution. From a statistical analysis of relaxed structures, we give 
evidence for four types of Ti-atom polar displacements: along the (111), (HO), or (100) directions 
of the cubic unit cell, or almost cancelled. The type of a Ti displacement is entirely determined by 
the Ti/Zr distribution in the adjacent unit cells. The underlying mechanism involves local strain 
effects that ensue from the difference in size between the Ti*+ and Zr cations. These results 
shed light on the structural mechanisms that lead to disordered Ti displacements in BaTii-^ZrajOa 
relaxors, and probably in other BaTiOa-based relaxors with homovalent substitution. 
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The understanding of ferroelectric materials is a very 
active research area of constant interest over the past 
decades and with great relevance to both fundamental 
and application-related issuesii Among ferroelectrics, re- 
laxors form a fascinating class of materials which exhibit 
complex structural and dynamical behaviors on multiple 
time- and length- scales^, with important technological 
applications such as capacitors and piezoelectric devices*^ 

Relaxors are characterized by a broad and frequency- 
dependent maximum of the dielectric permittivity as 
a function of temperature, in strong contrast with the 
sharp and frequency-independent anomalies observed in 
classic ferroelectrics. As the temperature decreases, most 
relaxors undergo no structural phase transition but de- 
velop a nanoscale structure that consists in polar atomic 
displacements with short-range correlations (so-called 
polar nanoregions). Chemical substitution was shown 
to be essential to obtain the relaxor properties, for it 
introduces spatial fluctuations of bonds and charges<^ 
The charge fluctuations, which occur in case of a het- 
erovalent substitution {e.g. Mg^^/Nb^"*" in the canonical 
relaxor PbMgi/3Nb2/303), were long considered as the 
main driving force to determine the polar nanoregions 
Nevertheless, the relaxor behavior is also observed in sev- 
eral lead-free BaTiOa-based solid solutions with a ho- 
movalent substitution, as exemplified by BaTii_j,Zrj;03 
(0.25 < X < 0.50).— Such systems are of great funda- 
mental interest, as they offer an unique opportunity to 
isolate and study the role of lattice deformations in the 
appearance of the relaxor behavior. 

In BaTii_j;Zr^03, the substituted Zr'*+ and Ti4+ 
cations exhibit very different sizes (r-ri4+ = 0.605 A, 
^Zr'i+ — 0-72 A) J Such a size mismatch was shown to 
induce distortions of the oxygen atoms network in the 
classic-ferroelectric perovskite PbZri_a;Tia;03, as well as 
a reorientation of the Pb polar displacements depend- 



ing on the Zr concentration!^ In BaTii_2;Zr2,03 relax- 
ors, the Zr substitution has to be thought as impeding 
the Ti polar displacement correlations that exist in the 
classic ferroelectric BaTi03, thus leading to the forma- 
tion of polar nanoregions. However, the exact structural 
mechanisms involved have not yet been conclusively elu- 
cidated in BaTii_2;Zra;03, nor in other BaTi03-based 
relaxors. "'^P"'^^ Earlier experimental studies of the local 
structure in BaTii_2;Zrj;03 showed that the Og octahe- 
dral cages around a Zr atom have a larger size (~4.20 A) 
than those around a Ti atom (~ 4.05 A) ^^'^^ , as expected 
from the different Ti*"''/Zr*+ cation sizes. In order to 
know if a Zr atom further affects its surrounding struc- 
ture, including the position of its Ti neighbors, one has to 
analyze the oxygen cage distortions and the Ti displace- 
ment changes as a function of the local Ti/Zr distribu- 
tion. This information cannot be found from experiments 
since the atomic positions are averaged out on all the Zr 
and/or Ti octahedral sites, even when using local probes 
such as x-ray absorption spectroscopy (XAS) or pair dis- 
tribution functions (PDFs). On the other hand, supercell 
first-principles calculations offer a promising (theoretical) 
route to study the local interplay between the Zr- and Ti- 
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To model Ti/Zr substitution effects in the 
BaTii_2;Zra;03 relaxors, we used three 135-atom 
supercells representative of BaTio.74Zro.26O3- The 
end-member compounds BaTi03 and BaZr03 were 
also considered as references to validate the calculation 
method. The initial atomic positions correspond to 
an assemblage of 3 x 3 x 3 perovskite unit cells in 
their ideal cubic form (space group Pm3m), with a 
lattice constant determined as the cube root of the 
experimental cell volume at 300 K (4.054, 4.006, and 
4.192 A for BaTio.74Zro.26O3, BaTi03, and BaZr03, 
respectivelyi^"— ) . For BaTio.74Zro.26O3, three pro- 
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FIG. 1: (Color online) Arrangements of the 20 Ti- and 7 
Zr-atoms in the BaTio.74Zro.26O3 supercells used for our cal- 
culations. Only one of the 27 perovskite unit cells is fully 
represented to facilitate the visualization of the supercells. 



totype configurations of the 20 Ti- and 7 Zr-atoms 
were considered. In the first supercell (Fig. [T^) the 
average number of Zr neighbors of a Ti atom is 1.5 
out of the 6 nearest octahedral B-sites, which is close 
to the value of 1.56 expected in the case of randomly 
distributed Zr atoms. The second (Fig. [T|d) and third 
(Fig. [T];) supercells were built in such a way that the 
average number of Zr neighbors of a Ti atom is 1.7 
(enhanced probability for Zr-Ti pairs of neighbors) and 
1.2 (segregation of the Ti and Zr atoms), respectively. 
Various geometrical arrangements of the Zr atoms were 
considered: triangular arrangements parallel to the (111) 
planes (supercell "a"), lines along the [100] directions 
(supercell "b"), square arrangements parallel to the 
(100) planes (supercell "c"), etc. 

Total energies calculations were carried out using 
the density functional theory as implemented in the 
Vienna ab-initio simulation package. Projected aug- 
mented plane wavesiS, (PAWs) with the Perdew-Wang 
exchange-correlation potential were adopted. The va- 
lence state of each element had been defined previously 
in the provided PAW potentials and the plane- wave cut- 
off was 400 eV. Numerical integrations in the Brillouin 
zone were performed by means of the Hermite-Gaussian 
method with N=l and smearing parameter of cr=0.1 eV. 
Brillouin-zone sampling using a 4 x 4 x 4 fc-point grid 
was found necessary for differences of total energies of 
our supercells to converge within 10""^ cV. Structural op- 
timizations were carried out under the conditions that 
all residual forces should be smaller than 0.01 eV/A. 
The BaTio.74Zro.26O3 supercell edges were constrained to 
form a cube with an edge-length of 3 x 4.054 = 12.162 A, 
in order to match the crystallographic structure deter- 
mined experimentally.— For BaTiOa and BaZrOa, only 
the supercell volume was fixed to the experimental value, 
the supercell-edge lengths and angles being free to vary. 

In the relaxed supercells of BaTiOa and BaZrOa, the 
atomic positions match within 1.5% those determined by 
neutron diffraction at low temperatures J^i^ We did not 
find any distortion with respect to the crystallographic 
structures, which is in agreement with the experimental 
studies of the local structure i^^i^^ This allows us to con- 
sider that the calculation method chosen is well suited for 
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FIG. 2: (Color online) Experimental pair distribution func- 
tion of BaTio.75Zro.25O3 at 300 K (open circles)^'*, compared 
with the PDF calculated from the three relaxed supercells of 
BaTio.74Zra.26O3 (line). 



studying these two compounds and, by extension, their 
sohd solution. 

Concerning BaTio.74Zro.26O3, the energies of the re- 
laxed structures are similar for the three supercells drawn 
in Fig. m within 0.31 eV (2.3 meV/atom). This en- 
ergy difference is small enough to consider that each of 
the "a" , "b" , and "c" supercells potentially represents a 
small region of the real material. We find that the in- 
teratomic distance distributions calculated for the three 
BaTio.74Zro.26O3 supercells are very similar, which sug- 
gests the existence of a single structural response to the 
Ti/Zr substitution, applicable for any local distribution 
of the Ti and Zr atoms. The PDF calculated from the 
summed interatomic distance distributions is compared 
in Fig. [5] with the experimentally determined neutron- 
PDF of BaTio.75Zro.25O3 at 300 K (Ref. [13).— An ex- 
cellent agreement is obtained between the two, especially 
in the r-range [1.5 - 2.5 A] that corresponds to the Ti-0 
and Zr-0 distances within oxygen octahedra, in which 
the most important deviations with respect to the crys- 
tallographic structure were observed^i^ Hence, the re- 
laxed supercells capture the essential features of the local 
structure in the BaTio.74Zro.26O3 relaxor, including the 
response of the Ti-Oe and Zr-Oe octahedra structures to 
the Ti/Zr substitution. 

We first determine the deformations undergone by the 
oxygen octahedra's network due to the Zr/Ti substitu- 
tion. Defining G as the center of mass of an octahe- 
dron, the octahedral distortions can be classified into two 
categories: angular distortions, which result in a tilt of 
the G-0 segments with respect to the Cartesian axes of 
the supercells, and radial distortions, which result in a 
change of the G-0 distances. The tilt angles of the G-0 
segments are found to be small, within the range [0 — 1 °] . 
They are not uniform in a given octahedron, and thus 
unambiguously correspond to angular distortions as op- 
posed to global rotations of the octahedra. No relation 
is found between the tilt angles and the local chemical 
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FIG. 3: (Color online) Distribution of the distance be- 
tween the center of mass of oxygens (denoted G) and 
the oxygens calculated from the three relaxed supercells of 
BaTio.74Zro.26O3: (a) for the Zr-Oe octahedra, and (b) for 
the Ti-Oe octahedra. The G-O segments that are directed 
toward Zr and Ti neighbors are represented in light magenta 
and dark green, respectively. The arrows mark the mean G-0 
distances calculated for the Zr-Oe and Ti-Oe octahedra. 



distribution. On the other hand, the G-0 distances ex- 
hibit a marked dependence on both the type (Ti-Oe or 
Zr-Oe) of the considered octahedron and the local Zr/Ti 
distribution on its neighboring sites. From the calculated 
G-0 distance distributions represented in Fig. [31 one can 
indeed observe that the average G-0 distance is larger in 
the Zr-Oe octahedra (2.092 A) than in the Ti-Oe octa- 
hedra (2.005 A). This result, which is expected from the 
different Ti*+/Zr'*+ cation sizes, is in quantitative agree- 
ment with the previous experimental studies of the local 
structureJ^i^ Besides, the large width of the peaks ob- 
served on the G-0 distance distributions suggests that 
the shape of the octahedral cages fluctuates from one 
site to the other. In order to understand whether the 
local Zr/Ti distribution determines or partly determines 
this disorder, we considered separately the G-0 distances 
that form a G-O-Ti chain (in dark green in Fig. [3]), and 
the G-0 distances that form a G-O-Zr chain (in light ma- 
genta in Fig. [3]). We find that the shortest G-0 distances 
are observed when the neighboring octahedral site is oc- 
cupied by a Zr atom, while the longest ones occur when 
the neighboring octahedral site is occupied by a Ti atom. 
The octahedra are thus compressed in the direction of Zr 
neighbors and expanded in the direction of Ti neighbors, 
which provides a microscopic demonstration of chemical 
pressure to which we will come back later. 

It can be expected that the aforementioned distortions 
of the oxygen atom's network modify the bonding of the 
cations in their octahedral and cubo-octahedral cages, 
and hence their polar displacement. In order to inves- 




FIG. 4; (Color online) Ti displacement vectors calculated in 
the three relaxed supercells of BaTio.74Zro.26O3, represented 
with the absolute value of their coordinates. Four types of dis- 
placement with distinct symmetries are observed: a) toward 
a face of the surrounding octahedron (full circles), b) toward 
its edge (open triangles), c) toward its corner (full triangles), 
and d) almost canceled (open circles). For comparison, the 
Ti displacement vector calculated in the classical ferroelectric 
BaTiOs is represented by a star. The symmetry of a Ti dis- 
placement is entirely determined by the distribution of Zr and 
Ti atoms in the adjacent unit cells, as illustrated in the inset. 
For the sake of clarity, only the Zr neighbors are represented 
for a given Ti06 unit. 



tigate this point we calculated the cation displacement 
vectors away from the center of mass of their surround- 
ing oxygen cages. The Ba and Zr atom displacements 
have a low average amplitude (0.06 A for Ba, 0.03 A for 
Zr), and their orientation is ill-defined. These charac- 
teristics indicate that the Ba and Zr atoms form mostly 
ionic bonds in the relaxor BaTio.74Zro.2603j as it is the 
case m the end-members BaTiOa and BaZrOs^i^ The 
average displacement amplitude calculated for the ferro- 
electrically active Ti atom in BaTio.74Zro.2eO3 is large 
(0.17 A) and nearly as high as the one calculated in the 
classical ferroelectric BaTiOa (0.18 A), which is in agree- 
ment with X-ray absorption near-edge structure studies 
of BaTii_^Zrj;03 at the Ti i^-edge^ However, the Ti 
displacement directions calculated for BaTio.74Zro.26O3 
and BaTiOa differ strongly. Figure S] gives a represen- 
tation of the calculated Ti displacement vectors, the x-, 
y-, and z-axes being parallel to the Cartesian axes of the 
supercell (i.e., the octahedra's axes before structural re- 
laxation). The Ti displacements in BaTiOs, which are 
all the same in amplitude and direction, are represented 
by the star. They have equal components along the x-, 
y-, and z-directions, and are thus directed along a (111) 
direction, toward a face of their oxygen octahedron. In 
the relaxed supercell of BaTio.74Zro.26O3, such Ti dis- 
placements toward an octahedron's face are still observed 
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(full circles), but other types of displacement also appear: 
along (110) toward an octahedron's edge (open triangles), 
along (100) toward an octahedron's corner (full trian- 
gles), and even nearly cancelled (open circles). Starting 
from a displacement toward an octahedron's face, the 
three other cases can be obtained by successively can- 
celling one displacement component along the x-, y-, or 
z-axes. By analyzing the Ti atom environments in the 
relaxed supercells of BaTio.74Zro.26O3, we establish that 
a Ti displacement component is cancelled as soon as a 
neighboring Zr atom lies on the corresponding octahe- 
dron's axis. This leads us to the important conclusion 
that in BaTio.74Zro.2603j the symmetry of a Ti-atom 
displacement is entirely determined by the distribution 
of the Ti and Zr atoms in the adjacent unit cells. Con- 
sidering the octahedral distortions caused by the Zr/Ti 
substitution, we can also conclude that a Ti atom is no 
longer displaced in the direction in which its octahedral 
cage is compressed. Hence, the strain induced by the 
Zr/Ti substitution leads to a reorganization of the cova- 
lent bonds between the Ti atoms and their oxygen neigh- 
bors in BaTio.74Zro.26O3: which results in a reorientation 
of the Ti polar displacements. 

It is instructive to set our results into the general con- 
text of the effect of strain in ferroelectric perovskites. 
More than thirty years ago. Samara et al. established 
that hydrostatic pressure (isotropic strain) reduces and 
even annihilates, at critical pressure, the correlations of 
polar cation displacements and thus ferroelectricityi22 In- 
dividual polar displacements are also sensitive to the 
pressure parameter. This is illustrated by XAS measure- 
ments on BaTiOa, which showed that the Ti polar dis- 
placements are reduced to zero due to the compression 
of the octahedra along their three axes<^ As a second 
illustration of the effect of strain, let us remind that com- 
pressive biaxial strain in perovskite thin films can greatly 
enhance the remanent polarization in the direction per- 



pendicular to the biaxial strain. Finally, we can cite the 
case of PbTiOs, where the Ti-Oe octahedra are elongated 
along one of their axes, due to the large displacement of 
Pb atoms. In this case, the local Ti displacements are 
directed in the direction of elongation, i.e. along a (100) 
direction rather than the preferred (111) ones.-^'"^^ It is 
interesting to note that the previously reported strain 
effects concern uniform distortions throughout the oxy- 
gen atoms' network, in opposition to our here reported 
random local strain effects induced by chemical disorder. 
Nevertheless, chemical pressure in BaTio.74Zro.26O3 acts 
on the individual polar displacements with similar rules 
to those observed for uniform lattice distortions. 

Our results provide insight into the structural mech- 
anisms that lead to disordered Ti displacements in 
BaTii_j;Zra;03 relaxors, and probably in other BaTiOs- 
bascd relaxors with a homovalent substitution. We show 
that the chemical substitution imposes a random dis- 
tribution of the Ti displacement symmetries among the 
Ti-06 octahedra, which impedes a perfect alignment of 
all the Ti displacements as it exists in the classic fer- 
roelectric BaTiOs at low temperaturesi^ Nevertheless, 
the local distribution of the substituted cations does not 
completely determine the Ti displacements: a displace- 
ment toward an octahedron's face can be directed in eight 
possible directions, four possible directions for a displace- 
ment toward an octahedron's edge, and two possible di- 
rections for a displacement toward an octahedron's cor- 
ner. This degree of freedom allows partial correlations of 
the polar displacements, and thus the formation of polar 
nanoregions. We believe that the reported behavior of 
the Ti polar displacements and its simple link to the lo- 
cal chemical order will give valuable constraints on both 
phenomenological and structural models which will aim 
at describing polar correlations in BaTi03-based relax- 
ors. 
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